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Abstract. The X-ray background intensity around Lick count galaxies and rich clusters of galaxies is investigated in 
three T^OS^Tenergy bands. It is found that the X-ray enhancements surrounding concentrations of galaxies exhibit 
significantly softer spectrum than the standard cluster emission and the average extragalactic background. The 
diffuse soft emission accompanying the galaxies is consistent with the thermal emission of the hot gas postulated 
first by the Cen & Ostriker hydrodynamic simulations. Our estimates of the gas temperature - although subject to 
large uncertainties - averaged over several Mpc scales are below 1 keV, which is substantially below the temperature 
of the intra-cluster gas, but consistent with temperatures predicted for the local intergalactic medium. It is pointed 
out that the planned ROSITA mission would be essential for our understanding of the diffuse thermal component 
of the background. 
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1. Intioduction 



20001 ). Its temperatures are between 10^ and 10^ K and 



The X-ray background (XRB) is mostly generated by dis- 
crete extragalactic sources (e.g. Lehmann et al. (2001) 



, and references therein). Among those sources, various 
$H classes of AGNs constitute a dominating part. Probably 
. 5^ , 5 to 10l% uf Lliu sufl XRD is piuduLud by huL gas in lIus- 
ters of galaxies. Aiuuud and beluw 1 kuV huL plasma in 
the Galaxy also contributes to the total background flux 
(Hasinger 1992| ). Apart from the source contribution, truly 
diffuse emission of extragalactic origin is also expected. 



the hottest and high density fragments surround mass con- 
centrations as groups and clusters of galaxies. Thus, the 
thermal emission by WHIM should be easiest to detect 
in the soft X-rays in the vicinity of rich clusters and high 



galaxy density areas. However, as pointed by Bryan & Voit 



(2001), compact sources produce most of the background 



Using hydrodynamic aimulationa Gen fe OatrikGi -( 1000 )- 



investi gated evolution of the primordial gao donoity and 
temperature. Baryons not condensed in stars and inter- 
stellar medium within galaxies, occupy the intergalactic 
space and are spread over a wide range of temperatures 



and observations put tight constraints on the level of truly 
diffuse XRB component. 

Detailed analysis of the diffuse emission produced by 
the intergalactic medium has been presented by Bryan & 
Voit (2001) and Croft et al. (2001). Using hydrodynamic 



and de n sities. This question is d i scusse d in d etail by Dave 
et al I (|2000|) , [Bryan fc VoiTI (|200l| ) and [Croft et al 



( 2001 ). According to all the simulations, the hottest phase 
is located in the high mass concentrations of clusters of 
galaxies and it is responsible for the cluster X-ray emis- 
sion. Tiplativply rn1d phase with lemppratnre below 1 fl^ K 
revealslits presence by the I.^fman alpha forest The frac- 



tion which still escapes detection, described as Warm-Hot 
Intergalactic Medium (WHIM), comprises of 30 - 40% 



simulations they generated maps of the X-ray sky which 
included contribution from the WHIM component. These 
maps have been used to calculate the model autocorrela- 
tion function of the XRB at small angular scales (below 
~ 10') as well as the cross-correlation with the galaxy 
distribution. 

In the present paper the distribution of the XRB in the 
vicinity of galaxy concentrations is carefully investigated 
from the observational point of view. Lick counts (Shane 
fc WirtanciT] [1967[, hereafter S'VVfl) an d Abell clusters of 
galaxies (Abell 195S , Abell et al. [ 1989[ ) have been used to 
select areas of the high mass concentrations. Since the pre- 



of all baryons in the present-day universe (Dave et al 



Send offprint requests to: A. M. Soltan 



^ The Lick galaxy counts in 10' pixels have been kindly pro- 
vided to us in the electronic form by Dr. M. Kurtz. 
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dieted WHIM temperatures are substantially lower than 
those of the intra-cluster gas, and the WHIM emission is 
softer that the average extragalactic XRB, the objective 
of our analysis is to look for the systematic variations of 
the XRB spectral slope as a function of the distance from 
peaks of the galaxy distribution. The expected WHIM 
contribution to the XRB is small in comparison to the 
total background flux. To increase signal-to-noise ratio we 
have measured the excess XRB flux around galaxy con- 
centration using the cross-correlation technique. The X- 
ray data contained in the ROSAT All-Sky Survey (RASS) 
have bee n used. For a comprehensiv e description of the 
RA SS see jSnowden fc Schmitt | ( |l990| ) , |Voges ct al. | ( |l999|) 



and [Snowden ct al. | ( |1995|) . Due to limited angular res- 
olution of the RASS we have con centrated on the la rger 
scales than those investigated by Croft et al. ( 2001 ). In 
the next section the basic concept of our investigation is 
presented and in Sect. ^ results of calculations are given. 
We end our investigation with the short conclusions in the 
Sect. HI 

2. The Method 

It is well established that the XRB exhibits fluctuations 



over a wide range of angular scales (e.g. Sliwa et al. ( 2001 ) 



and references therein). It is expected that the XRB vari- 
ations result from the nonuniform distribution of sources 
which generate the XRB (e.g. Soltan et al. | ( |1999| )). One 
should note, however, that intrinsically large, extended 
sources also produce fluctuations at the corresponding an- 
gular scales. Such fluctuations associated with rich clus- 
ters of galaxies have been detected by Soltan et al. ( 1996| ). 



They calculated the cross-correlation function (CCF) be- 
tween the distribution of the Abell clusters and the RASS 
maps. The region apparently free from the strong galac- 
tic emission has been used. The positive amplitude of the 
CCF extended up to ~ 7° and ^ 3° for nearby and more 
distant clusters, respectively. The CCF signal indicated 
that clusters are surrounded by halos of the enhanced X- 
ray emission. They estimated the average size of the halo 
at approximately 20Mpc. The amplitude of the X-ray en- 
hancements above the average XRB level is extremely low. 
It amounts to about 1 % of the total XRB surface bright- 
ness. Such weak signal could not be detected using an 
individual cluster observation but it was identified in the 
statistical sense in the cluster population using the cross- 
correlation analysis. 

Two explanations have been proposed in the original 
paper: discrete and diffuse. In the discrete model the ex- 
cess emission is produced by a large number of sources. 
The X-ray enhancements result from the fact that the 
spatial density of galaxies around clusters is larger. This 
increase d density extends up to a t least lOMpc from the 
clusters ( Seldner fc Peebles 1977 ). Such concentrations of 
galaxies naturally generate excess X-ray volume emissiv- 
ity, associated with the X-ray emission of normal galax- 
ies and weak AGN sources populating those galaxies. In 
the diffuse model it was ad hoc assumed that the X-ray 



emission was produced in the extended hot clouds of gas 
surrounding the clusters. 

These two models demonstrat e significantly d if ferent 
spe ctral proper t ies. A s shown by Cen fc Ostriker (1999) 



and pave et al. ( 200C ) the WHIM which would be a source 
of the X-ray enhancements contributes mostly at and be- 
low 1 keV, while the normal galaxy and AGN emission 
is substantially harder and is expected to have spectral 
shape not distinctly different from the integrated XRB 
flux. The hot intra-cluster gas with kT typically between 
2 and lOkcV also produces a harder spectrum than the 
WHIM. 

The CCF of the RASS and the SW counts has been de- 



termined in Soltan et al. (1997). The clear correlation sig- 
nal was interpreted as a result of the clustering of galaxies, 
but some indication of the extended emission associated 
with galaxies has been reported there. 

The excess emission Sp{9) at angular distance 6 from 
the randomly chosen galaxy or cluster is related to the am- 
plitude of the corresponding CCF, w{6), in the following 
way: 



dp{e) = {p)w{0), 



(1) 



where (p) is the average XRB signal. The details of the 
CCF calculations are given in ISohan et al.] ( |l996| ). Here 
we recall the main points. The final RASS data are binned 
into 12' X 12' pixels. We have binned the galaxy/cluster 
samples in the same way and the CCF was estimated using 
the formula: 



w{e) 



m{6) 



1, 



pn 



(2) 



where pi and nj denote the X-ray count rate in the i-th 
pixel and number of galaxies/clusters in the j-th pixel, re- 
spectively; m{9) is the number of pixel pairs with the sepa- 
ration equal to and the sum extends over all such pairs; p 
and n are the mean values of the XRB and galaxy/cluster 
distributions. 

The objective of the present calculations is to deter- 
mine the X-ray excess flux, Sp{9), as a function of the pho- 
ton energy. Unfortunately, the energy range of the ROSAT 
PSPC useful in the present investigation covers a rather 
narrow energy range, roughly between 0.7 and 2.0 keV ^ 
In the calculations we have used three standard ROSAT 
energy bands: R5, R6 and R7 with peak response at 0.8, 



1.1 and 1.5 keV, respectively (see Snowden ct al. 1994 for 
details). Since the PSPC energy resolution is relatively 
poor, those bands overlap what additionally reduces the 
sensitivity of our analysis. 

3. Results 

In the present calculations we have used the data at high 
galactic latitude in the northern hemisphere (6 > 40°), 



^ The soft band belo w 0.28 keV is highly contaminated by 
local galactic emission ( ^nowden et al. 1995). 
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where the extragalactic signal seems to be least contam 
inated by the galactic emission: 70° < I < 250° (3oltan 
et al. 1996). Although in this region some local contri- 



bution is still present, the CCFs between the RASS maps 
and optical samples are determined with relatively high 
accuracy. 

Since the present calculations are based on the an- 
gular cross-correlations, the effect depends on the galaxy 
distances. The redshift distribution of SW galaxies is rela- 
tively wide with a median value of 0.068 and 80 % of galax- 
ies have redshifts between 0.045 and 0.17. These estimates 
are based on extrapolations of the redshift distribution of 



the early SDSS data (|Dodelson ct al. ||200l| ) 



The AbeU clusters ( |AbeU ||195^ , |Abell et al.]|1989D have 



been divided into two distance groups. Sample I contains 
clusters with Distance Class 1 through 4 and all Richness 
Classes, and sample II contains DC = 5 and RC > 1. 
In the investigated area samples I and II consist of 57 
and 210 clusters, respectively. In our analysis both sam- 
ples give consistent results for all three energy bands (see 
below). However, due to the small number of clusters in 
the sample I, uncertainties of our measurements for that 
sample are too large to draw restrictive conclusions on the 
spectral properties of the X-ray halos. For this reason, the 
subsequent analysis is limited to the sample II. Using the 
Abell et al."| (|1989D data it is estimated that 80 % of DC 



= 5 clusters have redshifts between 0.07 and 0.17. 

The SW counts and Abell clusters cover compara- 
ble redshift ranges and the distributions of objects in 
both data sets are strongly correlated. To make the cross- 
correlation analysis of the SW and RASS data indepen- 
dent of the Abell clusters, all the calculations involving the 
SW galaxies were performed using a mask which removed 
the galaxies and X-ray flux related to clusters ^ 

The results of our calculations are shown in Figs. 1 
through 4 in the form of X-ray 'colour-colour' diagrams us- 
ing different combinations of bands R5, R6 and R7. To see 
effects of gas temperature variations with distance from 
galaxy concentrations, the excess fluxes have been calcu- 
lated using the corresponding CCF amplitudes at differ- 
ent separations. Data representing clusters are shown with 
squares while the data referring to SW galaxies - with 
pentagons. 

In Fig. the distribution of excess count rates in the 
R5 and R7 bands normalized to the cumulative count rate 
excess in all three bands is shown. The empty square rep- 
resenting the central cluster regions shows the emission de- 
termined from pixels containing the cluster center; the full 
symbol - represents the 'halo' region which is defined as 
an annulus with the inner and outer radius of 0?3 and 2?1. 
The solid and dotted ellipses denote 68% and 90% confi- 
dence areas. The SW data are shown for 4 separation bins: 
the largest symbol denotes signal for separations smaller 



Radius of the masked out area depends on the cluster dis- 
tance. For the nearest clusters (DC = 1) pixels within 48' from 
the cluster position are removed, while for DC > 5 only pixels 
centered on the cluster and 8 neighbouring pixels are not used. 



K 0.15 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 
R5/{R5 + R6+R7) 



Fig. 1. The X-ray 'colour-colour' diagram of the RASS 
data. Count rates in the energy bands R5 and R7 nor- 
malized to the summed rates in the bands R5, R6 and 
R7 are used. Open square indicates colours of the excess 
emission produced by clusters of galaxies; the full square 
- colours of the excess 'halo' emission integrated between 
0?3 and 2?1 from the clusters. The ellipses show 68% and 
90% confidence limits. Pentagons indicate colours of the 
excess emission correlated with the SW galaxies at 4 sep- 
aration bins: < 0?3, 0?3 - 0?7, 0?7 - 1?5, and 1?5 - 3?1, 
where the smaller symbol corresponds to the larger sep- 
aration; triangles show colours of AGNs in three redshift 
bins: 0.1-0.4 (large symbol), 0.4-1.0 (medium), 1.0-2.0 
(small), star - average RASS colours calculated from the 
count rates integrated over all the investigated area. Solid 
curve represents thermal Bremsstrahlung; labels indicate 
the temperature in keV. See text for details. 



than 0?3, and the next three pentagons of decreasing size 
represent three separation bins of 0?3 — 0?7, 0?7— 1?5, and 
1?5 — 3?1, respectively. 

To assess effectiveness and to inquire into potential 
systematic effects of our method, we have calculated also 
excess emission associated with the population of AGNs 



listed in the Veron catalogue (Vcron-Cetty 

pooo| ). 



Veron 



Thermal Bremsstrahlung spectra are given by the solid 
curve; labels denote temperatures in keV. In Fig. ^ the 
same data points are superimposed on the solid curve in- 
dicating the location of pure power law spectra; marks 
with labels give the photon indices. Potential effects of 
the Local Bubble and our own Galaxy halo emission on 
the integrated background flux are shown with two dashed 
lines. The lower line represents the sum of the power law 
spectrum with the photon index F = 2 and the thermal 
Bremsstrahlung with kT — 0.1 keV. The contribution of 
the thermal component increases along the line from at 
the left-hand end of the line to 0.1 of the counts in the 
R6 band at the right-hand end. The upper line represents 
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15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
R5/(R5 + R6 + R7) 



Fig. 2. All data points - same as Fig. |l|. Solid curve repre- 
sents colours of the pure power law models; labels denote 
the photon indices. Two dashed lines indicate loci of the 
two component models. The upper line: power law with 
index F = 1.4 with the addition (up to 10% in the R6 
band) of the thermal Bremsstrahlung with kT = 0.1 keV; 
the lower line - the same for F = 2. 



spectra with F = 1.4 and the thermal component varying 
in the same range. 

Uncertainty ellipses for the cluster sample have been 
delineated using the distributions of 1000 simulations. In 
each simulation, the count rates in randomly chosen RASS 
pixels were used to perform analogous computations as for 
the real data. In each simulated run the number of drawn 
pixels was equal to the number of clusters in the sample. 
Because the halo signal is substantially weaker than the 
'core' signal, our estimates of the halo temperatures are 
much less accurate. 

The error ellipses for the core and halo emission con- 
taining more than ^ 94 % of the simulation runs begin to 
overlap. This puts the lower limit for the significance of 
the temperature difference for these two cluster-centered 
regions. The statement is strengthened by the fact that the 
effect is present for any combination of the X-ray colours 
(see the remaining figures below). 

In the Figs. || and ^ the distribution of all the cate- 
gories of objects is shown using different combinations of 
the i?OS'^T bands. Although both SW and cluster points 
exhibit large uncertainties and scatter, all the data display 
a common feature. Namely, the emission parametrized by 
the temperature (or the spectral slope) correlated with the 
cluster and SW galaxies distribution becomes systemati- 
cally softer as the distance from the galaxy concentration 
increases. 

The results for the cluster sample can be summarized 
as follows: the cluster flux is consistent with the thermal 
emission with kT « 1.5 — 2keV (see below) and the halo 
flux - with kT « 0.5 — 0.8 keV. Although halo uncertain- 
ties are large and this temperature difference is barely 



significant, the general trend in all the 'colour-colour' di- 
agrams is highly suggestive. Similar temperature differ- 
ences between the core and halo excesses derived for all 
combinations of the energy bands demonstrate that the 
distribution of count rates in each band contributes to the 
resultant temperature effect around the cluster popula- 
tion. The excess associated with the cluster halos appears 
to be consistent with thermal emission generated by the 
WHIM. 

The excess emission correlated with the SW galaxies 
is still softer than that around clusters indicating lower 
gas temperature in qualitative agreement with simula- 
tions. Quantitative assessments of this effect are hindered 
by systematic effects (see below) as well as large statisti- 
cal errors. Although in the case of the SW points alone 
systematic colour variations with the increasing separa- 
tion are of low statistical significance, combined analysis 
of the cluster and SW samples shows that the effect is 
real. The error bars (not shown in figures for clarity) for 
three colours of the 'zero lag' SW data representing sepa- 
rations below 0?3 are in the range of 0.018 — 0.022. Thus, 
the SW signal is distinctly softer than the cluster 'core' 
emission. At larger separation the CCF amplitude suffers 
from substantially larger uncertainties and consequently 
colours of the '1?5 - 3?1' bin differ from the '< 0?3' bin by 
just over 1 a. One should note however, that the resultant 
CCF amplitude is a superposition of signals generated by 
the galaxy clustering and the extended emission surround- 
ing galaxies. This conclusion is supported by our earlier 
analysis of the SW - RASS correlation (see Boltan et al 



1997) 



Our estimates of the CCF between the RASS maps 
and the SW galaxies cannot be dir e ctly c ompared with the 
CCFs generated by |Croft et al.] (|200l|) for two reasons. 
Galaxies simulated by Croft et al. are more distant than 
the SW galaxies and, what is more important, those simu- 
lations are not devised to reproduce the discrete emission 
by galaxies, while in our measurements the CCF ampli- 
tude represents the total correlation signal. In the present 
investigation the IGM emission is revealed because it is 
substantially softer than the galactic emission. 

It is disturbing that our estimate of the average cluster 
temperature falls below 2keV. This value is significantly 
lower than the median temperature of w 5 keV observed in 



rich clusters (e.g. David et al. 1993, Horner et al. 1999) 



This discrepancy could be removed by a measured adjust- 
ment of the instrument calibration. However, modification 
of the effective area in bands R5, R6 and R7 which are re- 
quired to remove the cluster temperature inconsistency 
are of the order of 10 — 15 % what is larger than resid- 
ual uncertainties of the X-ray telescope/PSPC parame- 
ters (see below). In the recent study of the cluster X-ray 
emission distribution Bonamente et al. (2002) report de- 



tection of soft excesses in a large fraction of sources. The 
authors claim that in the soft ROS AT hand of ^ 0.25 keV 
detected fiuxes are on the average 9% larger than it is 
expected from extrapolation based on measurements at 
higher energies (above ~ 2keV). 
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0.2 0.25 0.3 0.35 0.4 0.45 0.5 
R5/(R5 + R6 + R7) 



Fig. 3. Same as Fig. |l| for the count rates in the bands R5 
and R6. 




0.35 0.4 0.45 0.5 0.55 0.6 

R6/(R5 + R6 + R7) 



Fig. 4. Same as Fig. |l| for the count rates in the bands R6 
and R7. 



To test if this phenomenon could explain our low 
temperature cluster measurement, we have calculated 
'colours' of the X-ray emission (defined by the R5, 
R6 and R7 bands) for the two temperature thermal 
Bremsstrahlung. Small amount of the low temperature 
emission was added to the main 'standard' component of 
5keV. It was found that ~ 0.7 keV emission normalized 



to 9 % in the low energy band used by Bonamente et al 



( 2002 affects the R5, R6 and R7 colours in such a way 
that they mimic the thermal emission with the tempera- 
ture of 2 — 2.5keV. Thus, this effect only partially explains 
our cluster temperatures, but the discrepancy is substan- 
tially reduced. 

A question of the 'soft excess' is complicated further 
by systematic uncertainties of the instrument calibration 
below the carbon edge (i.e. 0.28 keV). The comparison be- 
tween PSPC, Chandra and XMM soft sources shows that 



the PSPC gives about 20 — 30 % larger absolute fluxes 
at these energies. The agreement above 0.5 keV is much 
better (5 — 10%). Thus, the present investigation is prac- 
tically unaffected by uncertainties in the soft bands since 
the R5 - R7 bands have negligible responses below 0.5 keV 
([Snowden et al. 1 11994|) , but the [Bonamente et al. | (|2002D 



result could be subject to calibration errors. 

Systematic errors of the X-ray 'colours' in the RASS 
indicated by the cluster 'core' temperature are of the or- 
der of 0.05. Although we are unable to point out rigor- 
ously a single cause of this error, no potentially responsi- 
ble effect could alter the main conclusion of the present 
investigation, that the excess emission correlated with the 
galaxy distribution softens with the increasing size of the 
investigated area. Systematic errors in the 'colour - tem- 
perature' relationship introduce uncertainty of the cluster 
'halo' temperature at the level of 0.1 — 0.2 keV, and in 
the case of large separations in the SW galaxies - of about 
0.1 keV. Consequently, all our temperature determinations 
should be increased correspondingly. Thus, the cluster 
'core' temperature disagreement and statistical noise pre- 
cludes us from making a quantitative determination of 
the IGM temperature. But our investigation shows con- 
vincingly that concentrations of galaxies are surrounded 
by emission of significantly lower temperature than the 
average XRB and the cluster 'core' emission. 

One should note that the typical surface brightness of 
the IGM in the cluster halo is roughly two orders of magni- 
tude lower than the signal measured in pixels centered on 
clusters. Thus, a halo contribution to the amplitude of the 
CCF at zero lag is negligible. It implies that our estimate 
of the core temperature remains practically unaffected by 
the halo component. 

4. Summary and conclusions 

A perceptible fraction of the XRB fiuctuations is corre- 
lated with rich clusters of galaxies and with the overall 
galaxy distribution. The cross-correlations between the 
XRB and Abell clusters/SW galaxies extends up to several 
degrees. The large angular size enhancements of the XRB 
which surround regions of higher galaxy densities have dis- 
tinctly softer spectrum than the genuine cluster emission. 
It is highly suggestive that these enhancements are gener- 
ated by the WHIM. The halo emission is consistent with 
the thermal Bremsstrahlung with kT below 1 keV. Taking 
into account large uncertainties of the present estimates 
and a wide range of plasma temperatures predicted by Cen 
|fc Ostrikerl ( p^ ) and [Croft et al 



(2001), our results are 



in good agreement with theoretical models. Limited en- 
ergy range used in this investigation combined with the 
moderate energy resolution of the PSPC data prevent us 
from more detailed conclusions. 

This limitation applies also to our evaluations of the 
AGN spectral characteristics. Although the average X-ray 
colours of a population of AGNs are adequately mimicked 
by the thermal spectrum (Fig. |l|), simple power law spec- 
tral models also provide good fits (Fig. 0). Uncertainties 
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(not shown in figures) of the AGN colours are of the order 
of 0.009, 0.015, 0.050 for redshift bins 0.1 - 0.4, 0.4 - 1.0, 
1.0 — 2.0, respectively. Thus, all three samples do not ex- 
hibit significant colour differences. The majority of AGNs 
used in the present analysis are luminous X-ray sources 
and their average spectral properties are not identical with 
the average spectrum of weak sources which generate most 
of the XRB. Additionally, the integral soft XRB is con- 



taminated by thermal emission of our Galaxy (Hasinger 



1992). Importance of this effect in our calculations is il- 
lustrated in Fig. 1^ with dashed lines. 

Most if not all shortcomings of the present investiga- 
tion would be eliminated with the ROSITA mission which 
is being proposed. The ROSITA telescope mounted on 
board of the International Space Station would perform 
an all-sky survey within the energy range of 0.5 — lOkeV 
and a sensitivity about 100 times better than the RASS. 
These characteristics are particularly suitable for the anal- 
ysis of the low surface brightness features generated by the 
WHIM. Such investigation requires large, unlimited field 
of view. Also wide energy range extending well above the 
WHIM domain would help to isolate the WHIM thermal 
emission from the non-thermal components of the X-ray 
emission correlated with the galaxy distribution. High sen- 
sitivity would substantially increase efficiency of the cor- 
relation method by reducing the photon noise statistics. 

Finally, we would like to note that although the ex- 
tended clouds of hot gas around clusters in principle are 



source s of the Suny aev- Zel 'dovich (S-Z) ef fect (Sunyaev 
fc Zel ' dovich ( |198C| )), simple estimates by 3oltan et al 



( 1996 ) and Sliwa et al. (in preparation) indicate that for 
temperatures of the gas lower than 1 keV, the amplitude 
of the S-Z effect generated by halos is below the detection 
threshold of the COBE DMR measurements. 
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